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We propose a new approach to describe baryonic structure in terms of an effective chiral La- 
grangian. The state vector of a baryon is defined on the light front of general position (Ox = 0, 
i where (0 is an arbitrary light-like four vector. It is then decomposed in Fock components includ- 

ing an increasing number of pions. The maximal number of particles in the state vector is mapped 
i out to the order of decomposition of the chiral effective Lagrangian to have a consistent calcu- 

lation of both the state vector and the effective Lagrangian. An adequate Fock sector dependent 
renormalization scheme is used in order to restrict all contributions within the truncated Fock 
space. To illustrate our formalism, we calculate the anomalous magnetic moment of a fermion in 
the Yukawa model in the three-body truncation. We present perspectives opened by the use of a 
lO new regularization scheme based on the properties of fields as distributions acting on specific test 
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1. Introduction 

The calculation of baryon properties within the framework of chiral perturbation theory is the 
subject of active theoretical developments. Since the nucleon mass is not zero in the chiral limit, all 
momentum scales are involved in the calculation of baryon properties (like masses or electro-weak 
observables) beyond tree level. This is at variance with the meson sector for which a meaningfull 
power expansion of any physical amplitude can be done. While there is not much freedom, thanks 
to chiral symmetry, for the construction of the effective Lagrangian in Chiral Perturbation Theory 
(CPT), J??cpt, in terms of the pion field - or more precisely in terms of the U field defined by 
U = e n ' % ^ % where f K is the pion decay constant - one should settle an appropriate approximation 
scheme in order to calculate baryon properties. Up to now, two main strategies have been adopted. 
The first one is to force the bare (and hence the physical) nucleon mass to be infinite, in Heavy 
Baryon Chiral Perturbation Theory [1]. In this case, by construction, an expansion in characteristic 
momenta can de developped. The second one is to use a specific regularization scheme [2] in order 
to separate contributions which exhibit a meaningfull power expansion, and hide the other parts 
in appropriate counterterms. In both cases howether, the explicit calculation of baryon properties 
relies on an extra approximation in the sense that physical amplitudes are further calculated by 
expanding J^cpt in a finite number of pion field. 

Moreover, it has recently been realized that the contribution of pion-nucleon resonances, like 
the A and Roper resonances, may play an important role in the understanding of the nucleon prop- 
erties at low energies [3]. These resonances are just added "by hand" in the chiral effective La- 
grangian. This is also the case for the most important 2n resonances, like the a and p resonances. 

We would like to propose in the following a new formulation in order to describe baryon 
properties in a systematic way. Since in the chiral limit the pion mass is zero, any calculation of 
kN systems demands a relativistic framework to get, for instance, the right analytical properties of 
the physical amplitudes. The calculation of bound state systems, like a physical nucleon composed 
of a bare nucleon coupled to many pions, relies also on a non-perturbative eigenvalue equation. 
While the mass of the system can be determined in leading order from the iteration of the 7iN self- 
energy calculated in first order perturbation theory, as indicated in Fig.l.a, this is in general not 
possible in particular for nN irreducible contributions as shown on Fig.l.b. 




(a) (b) 

Figure 1: Iteration of the self-energy contribution in first order perturbation theory (a); irreducible contribu- 
tion to the bound state equation (b). 



2. Light-front chiral effective field theory 

The general framework to deal with these requirements is Light Front Dynamics (LFD). Rel- 
ativistic description admits some freedom in choosing the space-like hyper-surface on which the 
state vector is defined [4]. In the standard version of LFD, the state vector is defined on the plane 
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t + z = 0, invariant with respect to Lorentz boosts along the z axis. Since the vacuum state in LFD 
coincides with the free vacuum, one can construct any physical system in terms of combinations 
of free fields, i.e. the state vector is decomposed in a series of Fock sectors with an increasing 
number of constituents. Note that the triviality of the vacuum in LFD does not prevent from non- 
perturbative zero-mode contributions to field operators [5]. 

This decomposition of the state vector in a finite number of Fock components implies to con- 
sider an effective Lagrangian which enables all possible elementary couplings between the pion 
and the nucleon fields. This is indeed easy to achieve in chiral perturbation theory since each 
derivative of the U field involves one derivative of the pion field. In the chiral limit, the chiral ef- 
fective Lagrangian of order p involves p derivatives and at least p pion fields. In order to calculate 
the state vector in the Af-body truncation, with one fermion and (N — 1) pions, one has to include 
contributions up to 2 (N — 1) pion fields in the effective Lagrangian, as shown on Fig. 2. We thus 
should calculate the state vector in the Af-body truncation with an effective Lagrangian, denoted by 
<Sf£f f , and given by 



ce>N _ ot,p=2(N-l) n n 



For a non zero pion mass, one should extend this mapping by taking nr\ of order p = 2, as it is 
done in CPT. 



(N — 1) bosons 




(N — 1) bosons 



Figure 2: General vertex including a maximum of (N — 1) pion fields in the initial and final states. 



While the effective Lagrangian in Light-Front Chiral Effective Field Theory (LF# EFT) can 
be mapped out to the CPT Lagrangian of order p, the calculation of the state vector does not rely 
on any momentum decomposition. It relies only on an expansion in the number of pions in flight at 
a given light-front time. In other words, it relies on an expansion in the fluctuation time, xy, of such 
contribution. From general arguments, the more particles we have at a given light-front time, the 
smaller the fluctuation time is. At low energies, when all processes have characteristic interaction 
times larger than iy, this expansion should be meaningfull. 

It is interesting to illustrate the general features of LF# EFT calculations. At order N = 2, 
we already have to deal with irreducible contributions as shown on Fig.l.b. The calculation at 
order N = 3 explicitly incorporates contributions coming from nit interactions, as well as all low 
energy nN resonances, like the A or Roper resonances. Indeed, in the \nnN > Fock sector, the nN 
state can couple to both /=T = 3/2as well as/ = T = l/2 states. We can generate therefore 
all nN resonances in the intermediate state without the need to include them explicitly, provided 
the effective Lagrangian has the right dynamics to generate these resonances. This is the case, by 
construction, in CPT. 

To settle a general framework based on LFD, one has however to address three different prob- 
lems: i) one has to control in a systematic way the violation of rotational invariance; ii) one needs an 
adequate renormalization scheme consistent with the truncation of the Fock space; Hi) one should 
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use an appropriate systematic regularization scheme which preserves all symetries. We shall adress 
these three problems in the following. 

3. Covariant formulation of light-front dynamics 

In the Covariant Formulation of Light-Front Dynamics (CLFD) [6], the state vector is defined 
on the light-front plane of general orientation cox = 0, where CO is an arbitrary light-like four- vector 
co 2 =0. The state vector, 0© a (/?), of a bound system corresponds to definite values for the mass m, 
the four-momentum p, and the total angular momentum / with projection a onto the z axis in the 
rest frame, i.e., the state vector forms a representation of the Poincare group. The four-dimensional 
angular momentum operator, /, is represented as a sum of the free and interaction parts: 

J pv = Jpv ~t~ 'Jpv • (3.1) 
From the general transformation properties of both the state vector and the LF plane, we have: 

f p n lC(p)=L P A<0)<l>i O (p), (3-2) 

where 

L pv {co) = i[o> p ^-CO v ^ ) . (3.3) 

The equation (3.2) is called the angular condition. This equation does not contain the interaction 
Hamiltonian, once satisfies the Poincare group equations [6]. The construction of the wave 
functions of states with definite total angular momentum becomes therefore a purely kinematical 
problem. The dynamical dependence of the wave functions on the light-front plane orientation 
now turns into their explicit dependence on the four- vector CO. Such a separation, in a covariant 
way, of kinematical and dynamical transformations is a definite advantage of CLFD as compared 
to standard LFD on the plane t + z = 0. The eigenvalue equations for the Fock components can be 
obtained from the Poincare group equation [7] 

P<t>(p)=m 2 <t>(p). (3.4) 

We now decompose the state vector of a physical system in Fock sectors. Schematically, we have: 

(f){p) = \l) + \2) + ... + \n) + ... (3.5) 

Each term on the r.h.s. denotes a state with a fixed number of particles. It is proportional to the 
Fock component, or many-body wave function, (j) n . The spin structure of (f) n is very simple, since 
it is purely kinematical, but it should incorporate &>-dependent components in order to fulfill the 
angular condition (3.2). In the Yukawa model we have, for N = 2: 

= a\ u{k\)u{p) , (3.6) 



ri 2) = u(h 



m fb 

D2 



u(p) , (3.7) 



cop 

since no other independent spin structures can be constructed. In (3.6,3.7), we have used the vertex 
function r„ = (s n —M 2 )(j) n where s n = (k\ + . . .k n ) 2 . Here w's are free bispinors of mass m, a\, 
b\, and Z?2 are scalar functions determined by the dynamics. 
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4. Non-perturbative renormalization scheme in light-front dynamics 

In order to be able to make definite predictions for physical observables, one should also define 
a proper renormalization scheme. This should be done with care since the Fock decomposition of 
the state vector is truncated to a given order. Indeed, looking at Fig. 3 for the calculation of the 
fermion propagator in second order perturbation theory, one immediately realizes that the cancella- 
tion of divergences between the self-energy contribution (of order two in the Fock decomposition) 
and the fermion Mass Counterterm (MC) (of order one) involves two different Fock sectors [7]. 
This means that any MC and, more generally, any Bare Coupling Constant (BCC) should be asso- 

r r ^\ Sm 
+ — ^ 1 — + X 

Figure 3: Renormalization of the fermion propagator in second order perturbation theory. 



ciated with the number of particles present (or "in flight") in a given Fock sector. In other words, 
all MC's and BCC's must depend on the Fock sector under consideration. The original MC 8m 
and the fermion-boson BCC go should thus be extended to a whole series: 

<5m-+<5m« , g -> g$ (4.1) 

where / = 1,2, . . .N refers to the number of particles in "flight". A calculation of order N in- 
volves cW 1 ) . . . SmW and g^ . . .g^ . The quantities Sm^ and gff are calculated by solving the 
systems of equations for the vertex functions in the N = 1,2,3 .. . approximations successively. 
This procedure, which we call Fock Sector Dependent Renormalization, is a well defined, sys- 
tematic and non-perturbative scheme [7]. The MC is determined from Eq.(3.4), while the BCC 
g^ is determined by demanding that the co-independent part of the two-body vertex function T2 
at S2 = {k\ +&2) 2 = m 2 is proportional to the physical coupling constant g, b\(s2 = m 2 ) = g\fN[ 
where N\ is the normalization of the one-body Fock component (j>\ . 

5. The anomalous magnetic moment in the Yukawa model 

In order to address the calculation of a true non-perturbative system, we investigate the system 
composed of a fermion coupled to scalar bosons for the three-body, N = 3, Fock space truncation 
[8]. The system of equations one has to solve is given in Fig. 4, where the use of Fock sector 
dependent counterterms is shown. We use the Pauli-Villars (PV) regularization scheme, as detailed 
in [7]. The anomalous magnetic moment is calculated as a function of the boson PV mass, in the 
limit of infinite PV fermion mass. While the calculation for a boson-fermion coupling constant 
a = 0.2 shows nice convergence as a function of the PV mass, in the limit of large masses, the 
results for a = 0.5 are an indication that higher order Fock components may start to be sizeable. 

6. Perspectives 

It has been proposed recently [9, 10] to use the definition of quantum fields as operator valued 
distributions acting on test functions as a regularization/renormalization procedure. This construc- 
tion is particularly suited for light-front calculations since the test functions can be included from 
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Figure 4: System of equations for the vertex functions in the Yukawa model for the three-body Fock space 
truncation. Dashed lines correspond to scalar bosons. 
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Figure 5: The fermion anomalous magnetic moment as a function of the Pauli-Villars boson mass jl\ for the 
N = 3 Fock space truncation. We separate the contributions from the two-body (dashed line) and three-body 
(dotted line) vertex function to the total result (solid line). The lines are just drawn to guide the eyes. The 
results on the left correspond to a = 0.2 while the results on the right correspond to a = 0.5. The calculation 
are done with a fermion (m) and a boson (/i) mass of 1 GeV. 



the very begining in the definition of the Fock components, leading to a very transparent, and 
general, formulation [11]. Let us outline here the main steps of this formulation. The physical 
field (p(x) is defined in terms of the translation, T x , of the original distribution <I> acting on the test 
function p: 

<p(x) = TMP) = J d 4 yHy)p(x-y) - (6.1) 

We shall concentrate here on distributions singular in the UV domain. Any physical amplitude can 
be written in a schematic way (in one dimension for simplicity) as: 

s*= / dXT(X)f(X), (6.2) 
Jo 

where / is the Fourier transform of the test function, and T(X) is a distribution which may lead, 
without any regularization procedure, to a divergent integral in the UV domain when X —> <*>. In 
LFD, X is for instance proportional to k 2 , where k is the three-momentum of any of the con- 
stituents. The choice of the test function - of compact support and with all its derivatives equal 
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to zero at the boundary - should be done with care. In particular, one should make sure that the 
physical amplitudes are independent of the choice of the test function. This is achieved by using 
test functions which are partitions of unity, i.e. functions which are 1 everywhere except at the 
boundary where they should go to zero [10]. The choice of the boundary of the test function, H, 
should also respect scaling invariance which tells us that the limit X — > <x> can be reached in many 
different ways since in this limit r\ 2 X — > oo for an arbitrary scale T] 2 > 1. In order to do that, it 
is necessary to consider a boundary condition for which H depends on the kinematical variable X 
[10]. We can choose, for example: 

H(X) = ri 2 X a . (6.3) 

In the limit where a — > 1~, the maximum value of X, defined by H{X max ) = X max goes to infinity, 
and the test function goes to 1 in the whole kinematical domain. The extension of the distribution 
T(X) in the UV domain can thus be done easily using the general Lagrange formula given by 

m = ~I\ L 7 (1 [ xkf{xt) ] ' (64) 

for any integer k > 0. The physical amplitude writes in that case: 

tf= / dX^-^X k d* +l XT(X) / —(1-0* f(Xt)^ / dXT{X). (6.5) 

Jo k\ Ji t Jo 



in the limit where a — > 1~. This defines the extension t > {X) of the distribution T(X) in the UV 
domain. The value of k to be used depends on the degree of singularity of the initial distribution. 
The amplitude (6.5) is now completely finite. Note that we do not need the explicit form of the test 
function in the derivation of the extended distribution T > (X). We only rely on its mathematical 
properties and on its dynamical construction. 
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